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ABSTRACT 

The  response  of  semi- infinite  ice  floes  to  water  waves  is 
analyzed  for  relatively  deep  water.  If  the  floe  submergence  is 
neglected  it  is  found  that  a  progressive  wave  is  transmitted.  The 
stress  produced  by  this  transmitted  wave  is  determined  for  various 
floe  thicknesses  and  incident  wave  lengths.  When  the  submergence 
is  not  neglected  it  is  necessary  to  use  a  finite  difference  approach 
to  the  solution.  Such  a  solution  is  attempted  and  the  results  and 
accompanying  numerical  problems  are  considered  in  detail. 
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INTRODUCTION 


This  study  continues  the  effort  to  analyze  the  effect  of  water 
waves  on  ice  floes.  The  previous  study  (Ref.  1)  dealt  with  finite 
ice  floes  with  negligible  depth  of  submergence.  In  this  work  the 
effect  of  water  waves  on  semi- infinite  ice  floes  with  zero  and  finite 
submergence  depths  is  studied. 

The  first  portion  of  the  report  is  an  analysis  of  the  zero 
submergence  problem.  The  governing  equations  are  established 
and  a  closed-form  expression  for  the  stress  produced  by  the  trans¬ 
mitted  wave  is  obtained.  This  solution  depends  upon  the  roots  of  a 
quintic  equation,  which  involve  the  parameters  of  the  incident 
wave  as  well  as  ice  properties  and  floe  thickness.  This  stress  is 
calculated  for  various  floe  thicknesses  and  incident  wave  lengths 
assuming  constant  ice  floe  properties.  A  sample  calculation  is  in¬ 
cluded.  Reflection  and  transmission  coefficients  are  also  evaluated 
and  discussed. 

The  finite  submergence  case  is  then  analyzed.  Attempts  at  an 
analytical  solution  are  discussed  and  due  to  the  complexity  of  the 
equations  these  are  discarded. 

Finally,  a  finite  difference  approach  to  the  solution  is  used 
to  determine  the  stress  in  the  floe  near  the  leading  edge.  The 
limitations  of  computer  storage  severely  restrict  the  size  of  the 
flow  field  under  consideration  but  a  definite  indication  to  the  stress 
pattern  is  established. 
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SEMI-INFINITE  PLATE  NO  SUBMERGENCE 

Consiier  a  semi-infinite  elastic  plate  of  thickness  h  float¬ 
ing  upon  water  of  relatively  infinite  depth!  as  shown  in  Fig.  1. 

If  linearized  theory  is  used  then  the  free  surface  portion 
of  the  flow  field.  Region  1,  is  given  by  the  potential  function 

r  /T>  -ikx  .  „  ikx.  -ky  itrt 

(1)  $1  =  (Be  +  R  e  )  e  7  e 

with 

*ii  9M  i  9M 

(2)  TT  *  -By)  y  =  0  ’  ^1  =  '  i  SC)  y  =  0 
and  therefore 

(3)  «r2  =  g  k 

where 

=  local  elevation  of  profile  (ft.  ) 

| rj ^  |  =  incident  wave  amplitude  (ft.  ) 

Lj  =  incident  wave  length  (ft.  ) 

Tj  =  incident  wave  period  (sec.  ) 

<r  =  2ir/T1 

k  =  2w/L1 

U  =  hj  |g/«r 

R  =  reflection  coefficient 

B  =  i  U 
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The  equation  for  the  oscillation  of  a  plate  is 


(4) 


D 


Psh 


P 


where 


D 


h 

ps 

p 

v 

E 


displacement  of  the  plate 


h3  E 


12(1-  v  ) 


flexural  rigidity  of  the  plate 


plate  thickness 
mass  density  of  plate 

pressure  exerted  on  the  plate  by  the  fluid 
Poisson's  ratio  of  the  plate  =  i- 
modulus  of  elasticity  of  the  plate 


Now  if  $2  *s  *ke  potential  function  for  Region  2  and  if 
is  of  the  form 


(5) 


*2 


A 

n 


iX  x 
n 
e 


e 


X 

n 


y 


i<rt 

e 


then  the  pressure  p  as  shown  in  Ref.  1  is  given  by 


(6) 


P 


"  pf 


3*2  \ 

ST)  y  =  0  +  8  ^2 


where  is  the  fluid  density. 

If  the  plate  displacement  is  equal  to  the  water  displacement 
and  the  velocities  correspond,  then 


-4- 


dr\- 


9§, 


9t  ~  9y  j  y  =  0 


or 


^2  = 


i  aM 

cr  3y  /  y  =  0 


Thus,  the  plate  displacement  is 


(7) 


x\  x  -X  y  i<rt 

1  \  a  n  n7 

ri_  =  —  X  A  e  e  e 

'2  cr  n  n 


where  the  X  are  five  roots  of  the  characteristic  equation  (9),  which 
n 

is  obtained  by  substituting  (7)  into  (4) 


(8) 


i 

£T 


dx 


-  9a^(T  +  Pfg 


8$. 


•)yo  '  ‘ 


io"Pf  §2  =  0 


But 


8$_  \  iX  x  i<rt 

— a —  )  =  -  X  A  e  n  e 

8y  /  y  =  0  n  n 


therefore  the  characteristic  equation  is 


(9) 


X  +  H  (1-ak)  X  -  Hk  =  0 

n  '  '  n 


where 


H  = 


Pfg 


P„  h 


a  = 


and  k  =  o*  /g  for  the  incident  wave. 

The  roots  of  the  characteristic  equation  are  X^  .  With  a,  b,  c,  d 
all  greater  than  zero  the  X^  may  be  written  in  the  form 
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Xq  a  positive  real  number 
=  a  +  bi  X^  =  a  -  bi 

X-  =  -  c  +  di  \  A  =  -  c  -  di 

3  4 

but  since  the  potential  function  J  must  decay  with  increasing 

y,  X^  and  X ^  are  not  admissible.  The  potential  can  t>e 

written  as 

/  -iX  x-X  y  iX.x-X.y  -iX  x-\  yv  i<rt 

(10)  §2  =  (Te  °  °  +  Al  e  1  +A2e  2  2  )e 

(-iXQX-Xoy) 

Now  e  represents  a  progressive  wave  as  x 

approaches  infinity  and  thus,  the  term  T  is  related  to  the  trans- 

(iX.jX-X.jy)  (-iX.  x-X  y) 

mission  coefficient.  The  terms  e  and  e 

are  waves  which  decay  both  with  increasing  x  and  increasing  y  . 

Now  for  continuity  along  the  line  x  =  0  the  potential 

functions  and  particle  velocities  must  be  identical  so  that 

=  at  x  =  0 

—  at  x  —  0 

The  plate  has  a  free  end  at  x  =  0  so  that  the  moment  and  shear 
must  be  zero  at  that  point,  giving 

*2  a3 

d  d  n? 

(12)  - 1  -  0  at  X  =  0  ,  - 1  =  0  at  x  =  0 

dx  dx 3 


*1 

"33T 


-6- 


Now  (11)  and  (12)  give  four  boundary  conditions  which  are 
sufficient  to  determine  the  four  unknown  coefficients  R,  T, 
and  * 

If  the  subscript  (  )^  refers  to  the  real  part  of  a  coefficient 
and  (  )^  refers  to  the  imaginary  part  of  a  coefficient  (T^  is  not  to 
be  confused  with  the  period  of  the  incident  wave,  which  is  not 

used  explicitly  in  the  analysis),  then  after  considerable  algebraic 
manipulation  it  is  found  that 


T1 

= 

-  Zy  A  GU/J 

T2 

= 

2\  A  FU/J 

A11 

= 

IS  [*v3'v4*  Tl  + 

<V64>  T2] 

A12 

= 

'IS  [  <63-64>  T1 

1 

j 

O J 

I 

H 

A21 

= 

-i  [<WT1 

-(63*64)  T2 

A22 

= 

"  IS  [  (63*64>  T1 

*(v3  +  y4)  t2 

R1 

= 

T1  +  A1 J  +  A21 

R2 

T2  +  A12  +  A22  ‘  U 

where 


=  Au  +  i  A 


12 


=  A21  +  i  A 


22 


V 


6 

n 


etc.  and 


A  =  >3V4  +  6364 

F  =  A  (1  +y)  -  \y4  *  V3  ~  &3 


-7- 


G 


=  V*3+  V4‘61*4 


J  =  F2  +  G2 


In  order  to  find  the  relation  between  the  transmitted,  reflected 
and  input  waves  define  EE  =  RR  +  TT  where  the  bar  indi¬ 
cates  complex  conjugate,  then 


(14)  EE  =  R,2  +  3/  +  Tj2  +  T22 


Upon  substitution  of  the  quantities  T.  and  R^  from  (13)  it  is 
found  that 

"5>  -  i^[i)42-Mn(-i-i-y,*v(l 

♦YY44Y,Y34  Vj] 

Now  EE  is  a  positive  quantity,  being  a  sum  of  squares  of  real 

E  E 

numbers,  therefore  — exceeds  unity.  This  is  an  apparent 

U 

violation  of  the  equation  of  conservation  of  energy  which  apparently 

would  require  E  E  =  =  R  R  +  T  T  .  But  work  has  also  been 

done  in  permanently  changing  the  wave  length  from  =  2ir/k 

for  the  incident  waves  to  L  =  2  W  X.  for  the  transmitted  waves 

o  o 

and  (15)  expresses  this. 

Fig.  2  shows  the  variation  of 

_  (EE) 1/2  _  (R  R) 1/2  _  (T  T) 1/2 

t  -  - E -  *  P  -  - 0 -  ■  t  -  - u - 


8- 


FIGURE  2 

REFLECTION  AND  TRANSMISSION  COEFFICIENTS  VS.  INCIDENT  WAVE  LENGTH 


with  incident  wave  length  and  plate  thickness* 

The  stress  in  the  plate  reaches  a  maximum  value  near  the 
leading  edge  and  then  decays  to  a  slightly  lower  value  given  by  the 
amplitude  of  the  transmitted  wave.  The  peak  stress  is  analogous 
to  a  transient  and  the  stress  due  to  the  transmitted  wave  is  analogous 
to  a  steady  state  value  in  the  sense  that  it  is  independent  of  x  . 

The  transient  stress  is  found  from  the  expression 


<=u2  ♦  «=2*2  )  “ 


where 


E  hi 
s  o 

(1  -  v2) 


<T1  P1  *  w 


where 


2x 


E  hX 


Zy  <r(l-v  ) 


(■TlP2*  W 


-bx, 


=  -  Vj  cos  X^x  +  e  (y^  cos  ax  +  6^  sin  ax) 


=  -  y .  sin  X  x  +  e 


-bx  . 


sin  ax 


and  the  units  of  will  be  lb/ft  per  foot  of  incident  wave  amplitude. 

The  stress  due  to  the  transmitted  wave  S*  is  found  by  allowing 
x  to  approach  infinity  in  the  preceding  expressions 

2  \  1/2 


S* 


(C.2  *  C22) 


where 


E  hX 


2  cr  ( 1  -v  ) 


(-T.  cos  X  x  -  T,  sin  X  x) 
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T-,  cos  X  x) 
2  o 


Eh  \ 

8  o 


—  (T.  sin  \  x  — 

2.  '  1  o 


giving  finally 
S* 


2(r{l-v  ) 


9EhXoJYAU  ^ 
- -  J 


for  v  =  1/3. 

The  quantity  S*  for  various  incident  wave  lengths  and  plate 
thicknesses  is  shown  in  Fig.  3  . 

For  comparison  with  the  finite  submergence  case,  the  stress 
amplitude  for  an  incident  wave  length  of  10  feet  and  plate  thick¬ 
ness  of  5  feet  has  been  worked  out  in  detail  as  shown  in  Fig.  4  . 

As  an  illustrative  example,  suppose  it  is  desired  to  find 
the  stress  amplitude  for  a  15  foot  thick  floe  at  an  incident  wave 
length  of  100  feet  with  an  incident  wave  amplitude  of  6  inches. 

From  Fig.  3,  S*  for  h=  15,  =  100  is  112  lb/in^  therefore 

the  desired  value  of  S*  will  be 


S*  =  yj  112  =  56  lb/in2 


since  the  stress  amplitude  graph  is  normalized  to  a  1  foot  amplitude. 
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FIGURE  4 

STRESS  AMPLITUDE  IN  PLATE  (NO  SUBMERGENCE) 
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SEMI-INFINITE  PLATE  WITH  FINITE  SUBMERGENCE 


We  consider  a  semi-infinite  elastic  plate  submerged  in  water 
to  a  depth  "a 11  below  the  mean  water  line  of  a  train  of  incident 
harmonic  waves.  The  coordinate  y  is  measured  positive  down¬ 
ward  from  the  mean  water  line  and  the  coordinate  x  is  measured 
parallel  to  the  mean  water  line  from  the  leading  edge  of  the  plate 
positive  in  the  direction  of  incoming  wave  propagation  (i.  e.  positive 
into  the  plate). 

The  free  surface  portion  of  the  flow  field  must  obey  equation 
(2)  of  the  previous  section  of  this  report.  It  is  assumed  that  this 
submerged  edge  of  the  plate  perturbs  a  relatively  small  portion  of 
the  flow  field  and  therefore  that  as  the  distance  from  this  area 
increases!  the  general  simple  harmonic  character  of  the  incident, 
reflected  and  transmitted  waves  is  maintained. 

Ur  sell  (Ref.  2)  has  shown  that  if  a  normal  velocity  is  pre¬ 
scribed  on  part  of  a  vertical  plane  (i.e.  thin  boundary)  extending 
from  the  surface  then  the  two  dimensional  problem  of  determining 
the  motion  on  either  side  of  the  plane  as  well  as  beneath  it  can  be 
solved  for  the  deep  water  case.  The  motion  is  determined  by  a 
comparatively  simple  integral  equation.  At  the  lower  edge  of  the 
plane  the  velocity  becomes  infinite  due  to  a  mathematical  singu¬ 
larity  of  logarithmic  type.  Ur  sell  obtains  his  solution  by  assuming 
that  the  transmitted  waves  have  the  same  wave  length  as  the  incident 
waves. 
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FIGURE  5 

FINITE  DIFFERENCE  APPROXIMATIC 
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DN  TO  THE  FLOW  FIELD 
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An  attempt  was  made  to  modify  Ur  sell's  approach  to  include 
the  effect  of  the  plate  as  well  as  the  change  in  wave  length  of  trans¬ 
mitted  waves.  In  so  doing  a  very  complicated  set  of  simultaneous 
integral  equations  is  obtained.  The  kernels  of  these  equations 
involve  the  fifth  order  equation  for  the  plate  oscillation.  The  inte¬ 
grals  were  found  to  be  intractable  and  hence  the  analytical  approach 
was  discarded  in  favor  of  a  numerical  means  of  solution. 

Since  no  analytical  solution  was  found  for  the  plate  with 
finite  submergence  an  approximate  solution  was  attempted.  The 
coordinate  system  is  shown  in  Fig.  5. 

The  incident  waves  are  given  by  the  potential  function 

(16)  =  (  iUe”ikx  +  Re1**)  e*^  ei<rt 


The  submerged  edge  of  the  plate  is  assumed  to  perturb  a 
relatively  small  portion  of  the  flow  field,  therefore,  at  some 
station  AB  where  x  =  ,  the  potential  is 

(-ikx.  ikx.  \  -ky  i<rt 

iUe  1  +  Re  Aj  e  e 


Now 


d$ 

as  y  increases  -g-=-  vanishes,  thus,  for  sufficiently  large  y 


(18) 


8*BC 

~W~ 


y  =  y. 


=  o 


Again  for  sufficiently  large  x  =  x^  the  potential  which  is 
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the  asymptotic  form  of  given  in  equation  (10)  is  assumed  to 
exist  in  the  form 

-iX  x  -X  y  i<rt 

(19)  $CD  =  Te  °  e  °  e 


Along  the  boundary  DE  the  equation  of  motion  of  the  plate  must 
hold  and  this  is  given  by 


(20) 


+  H(l-ak) 


9x 


8$ 


DE 


8y 


,  ♦  Hk  (*de)  -  0 

y  =  a  '  7  y  =  a 


Along  the  submerged  edge  of  the  plate  EF 


(21) 


a*EF 
~~ 3x“ 


0 


since  this  edge  does  not  move  in  the  x  direction.  The  free 
surface  condition  along  FA  is  analagous  to  equation  (2) 


(22) 


$  = 


y  =  0 


Since  the  plate  has  a  free  end  at  the  point  x  =  0,  y  =  a, 


(23) 
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The  field  equation  is  Laplace's  Equation 


(24)  V2#  =  0 

which  holds  everywhere  within  A  B  C  D  E  F. 

There  are  now  sufficient  equations  to  determine  $  ,  R  and  T. 

The  size  of  the  grid  i  for  a  finite  difference  solution  of 
the  problem  is  unfortunately  not  arbitrary.  The  core  storage 
capacity  of  the  presently  available  digital  computer  equipment 
(IBM  7090)  limits  the  size  of  the  field  considerably. 

The  following  parameters  are  used  in  the  computation 
=  5  ft. 

=  1.44  x  108  lb/ft2 

=  1/3 

=  1 . 025  gm/cm^ 

=  0.9  Pf 

=  4.  5  feet 

=  %  =  1.125  ft. 

Equal  divisions  of  length  i  define  the  nodes  in  the  free- 
surface  portion  of  the  field.  There  are  10  nodes  in  the  y-direction 
and  at  most  8  in  the  x-direction.  The  field  under  the  plate  is  a 
semi-logarithmic  distribution  of  nodes  where  the  logarithmic  scale 
is  in  the  direction  of  increasing  x  and  the  y  divisions  are  I 
apart  as  in  the  free-surface  region.  There  are  16  nodes  in  the 
x-direction  and  6  in  the  y-direction. 
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The  equation  for  the  x-divisions  is 


<25>  *  *  rog-zTTiVzs  108  27  '  108  (27-m) 

where  0  <  m  <  1 6  in  unit  steps. 

In  order  to  conserve  computer  storage  and  yet  include  as 
much  of  the  perturbed  flow  field  as  possible,  the  free- surface 
portion  of  the  field  was  stepped  as  shown  in  Fig.  5  . 

Assuming  a  second  order  function  for  <|>  then  the  field 
equation  can  be  written  as 


v2# 


(26) 


S(«  +  yT  §1  +  *2  +  7T^+7T  *3  +  “7  *4 

-(=•?)  *•  ■  • 


where  a  and  y  and  are  defined  in  Fig.  6 

Along  the  free  surface  FA  the  slope  condition  (22)  permits 
elimination  of  the  ficticious  point  "1  M  giving,  since  a  =  y  =  i  , 


?2*FA  -  jp.  [»l+  2*4  +  *3-  H-2W)i0] 


=  0 


where  the  $(s  have  the  same  relations  to  as  shown  in  Fig.  6 
Along  the  lower  boundaries  AB,  BC,  CD  Eqn.  (24)  the 
field  equation  becomes  respectively, 


(28) 


^7  (*2  +  *3  +  *4  •  4*o) 
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FIGURE  6 

FINITE  DIFFERENCE  SCHEME  FOR  THE  FIELD  EQUATION 
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where  is  given  by  (17) 


(29)  V2§bc  =  a(Av>  *l*JZ  *i+  TT^*y)  *3  ‘  (^  *  71)  *°  '' 
(30>  v2iCD  -  *,  h**4  -{■Z*7z)io-  ° 


Along  the  boundary  DE  the  high  order  differential  equation 
is  expressed  in  finite  difference  form  using  the  scheme  shown  in 


Fig.  7 


Defining  ZN  =  -  x^  and  y'  =  y  -  4 f  then,  if 


3§  \ 
*7)  ^ 


y'  =  0 


and  VN  is  a  fourth  order  equation  in  x 


*n+6  -  Fn 


l  M-z"" 


The  constants  can  be  evaluated  by  curve  fitting  as  before.  It 

is  found  that 


hf  =:  vV-vV  +  v  V  -  v  V  +  V  V 
4  1  1  2  2  3  3  4  4  5  5 


where 


"2  “3  4  5 
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FIGURE  7 

DETAIL  FOR  EQUATION  (20)  AT  POINT  "1"  (SCHEMATIC) 


VZ  ■  ZjiVVV  z2mz3-z2) 

v3  =  z3(Z5-23)(Z4- Z3)(Z3-Z2) 

V4  *  24f25-  Z4l(24L-"Z-3rCZ'4-  z2) 

y5  =  Z5(Z5-  z4mz5'-ZjMZ5-  Z2) 


Now  using  the  field  equation  =  0  the  ficticious  points 

in  (31)  can  be  eliminated.  For  example,  at  point  "n" 


so  that 


(33) 


*N-I 


<XN-XN-l><XN+l-XN-l)  2<' 


(xn-xn-iHxn+i'Xn) 
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Finally,  equation  (20)  becomes 


(34)  24  M.  +  H(l-ak)  V  +  Hk  $  =0 

4  N  N 

Similar  equations  express  the  conditions  at  the  free  end  (23). 
Along  the  boundary  EF  using  (21) 

(35)  -L  (2§j  +  *2+  §4-  4i0)  =0 

The  system  illustrated  in  Fig.  5  consists  of  332  unknown 
values  of  §  where  4>^n  j  =  Re  (§)  and  <}>  ^  =  Im  ($),  plus  four 

additional  unknowns:  The  transmission  coefficient  T  =  Tj  +  i  T^ 
and  reflection  coefficient  R  =  +  iR^.  There  are  336  equations 

in  the  system;  332  statements  at  node  points  and  four  additional 
equations  for  the  free-end  conditions  at  x  =  0. 

The  value  of  the  stress  amplitude  is  found  from  the 

equation 

(36)  Sx  =  sts/  +  S2Z)l/Z 
where 


E  h 
8 


and  s  =  - 5“ 

2(1-  v  ) 
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Values  of  the  stress  are  shown  in  Fig.  8,  where  the 

corresponding  values  determined  for  the  same  plate  without  sub¬ 
mergence  are  also  shown. 

In  general  the  solution  of  the  finite  difference  formulation 
of  the  problem  of  non-zero  submergence  involves  the  following 
steps: 

1.  Establish  the  finite  difference  grid  and  idei  «,ify  the 
nodes. 

a.  Choose  i  . 

b.  Choose  the  semi-logarithmic  law  analagous  to 
equation  (25)  for  the  variation  of  x  beneath  the 
plate. 

2.  Set-up  the  linear  equations  at  each  node.  There  will 

be  two  equations  at  each  node  (except  at  the  plate  corner 
E).  One  equation  will  be  for  the  Real  part  of  §  and 
the  other  for  the  imaginary  part. 

a.  Use  Eqn.  (26)  at  each  interior  node. 

b.  Use  Eqn.  (27)  along  the  boundary  FA. 

c.  Use  Eqn.  (28)  along  the  boundary  AB. 

d.  Use  Eqn.  (29)  along  the  boundary  BC. 

e.  Use  Eqn.  (30)  along  the  boundary  CD. 

f.  Use  Eqn.  (34)  along  the  boundary  DE. 

g.  Use  Eqn.  (35)  along  the  boundary  EF. 

h.  At  point  E  state  equations  (23)  using  and 

from  Eqn.  (31). 
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FIGURE  8 

STRESS  AMPLITUDE  IN  PLATE 


Band  the  matrix  and  solve.  The  computer  routine 
used  for  this  report  was  "Linear  Equation  Solver  of 
Banded  Matrices",  SHARE  Distribution  Number  990 
CORR/1490  F4RWLE4F  from  the  SHARE  library. 


Using  the  solution  for  Re  §  and  Im  $  determine 

^7 

given  by  Eqn.  (36). 


-g—  from  Eqn.  (33),  then  finally  the  stress  will  be 
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ANALYSIS  OF  RESULTS  AND  CONCLUSIONS 


The  interaction  of  progressive  waves  with  a  semi- infinite 
ice  floe  has  been  investigated  wherein  the  ice  floe  is  considered 
equivalent  to  an  elastic  sheet.  In  this  investigation  two  approaches 
have  been  used;  neglect  of  the  finite  submergence  of  the  floe  and 
consideration  thereof.  In  the  case  of  the  former  approach,  a  mathe¬ 
matical  method  is  developed  wherein  a  closed  form  analytic  solu¬ 
tion  for  the  ice  floe  response  may  be  obtained.  The  results  of  the 
mathematical  analysis  show  that  a  portion  of  the  incident  wave 
energy  is  reflected  while  some  energy  is  transmitted  beneath  the 
ice  in  the  form  of  a  progressive  wave  whose  frequency  is  dependent 
on  both  the  characteristics  of  the  ice  and  the  incident  wave  length. 
The  partition  of  wave  energy  into  transmitted  and  reflected  energy 
is  represented  by  potential  function  coefficients  denoted  as  trans¬ 
mission  and  reflection  coefficients.  The  variation  of  these  coeffi¬ 
cients  as  a  function  of  incident  wave  length  for  a  5  foot  thick  floe 
is  shown  in  Figure  2. 

Further,  the  absolute  value  of  the  bending  stress  in  the  ice 
floe  due  to  the  transmitted  wave  is  shown  in  Figure  3  as  a  function 
of  floe  thickness  and  incident  wave  length. 

A  specific  example  showing  the  variation  of  the  bending  stress 
with  the  distance  from  the  free  edge  of  the  ice  floe  is  shown  in 
Figure  4.  The  parameters  chosen  for  the  calculation  were  an  ice 
thickness  of  5  feet  and  an  incident  wave  length  of  10  feet  (also  a 
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unity  amplitude  for  the  incident  wave  height).  When  observing  the 
stress  variations  depicted  in  Figures  3  and  4  certain  salient  points 
should  be  observed.  First  the  stresses  shown  represent  "absolute1' 
values.  However,  the  variation  of  the  ice  response  is  harmonic 
in  time.  Hence,  at  a  given  location  in  the  ice,  it  is  clear  that  the 
stress  values  shown  in  these  figures  alternate  between  the  positive 
value  shown  and  the  negative  of  this  value.  Secondly,  in  referring 
to  Figure  4  which  shows  the  stress  variation  as  a  function  of 
position  from  the  free  edge  of  the  floe,  it  is  important  to  realize 
that  the  "transient"  stress  values  (i.e.  those  values  close  to  the 
free  edge)  are  within  about  10  per  cent  of  the  "steady  state"  stress 
values  (at  positions  far  removed  from  the  free  edge).  Also  note 
that  the  transient  effect  exists  for  distances  greater  than  10  wave 
lengths  of  the  incident  wave  from  the  free  edge.  Although  only  one 
specific  example  has  been  calculated  for  the  transient  response, 
this  example  is  important  from  a  comparison  standpoint  with  the 
numerical  example  worked  for  the  same  floe  but  where  the  finite 
submergence  has  been  included. 

Turning  now  to  the  numerical  example,  some  comment  on  the 
choice  of  parameters  and  the  mathematics  behind  the  results  are  in 
order. 

It  is  clear  that  the  choice  of  an  input  wave  length  of  10  feet 
and  an  ice  thickness  of  5  feet  are  physically  unrealistic  (Ref.  3).  The 
necessity  for  the  choice  of  parameters  is  due  to  the  computational 
procedure  that  we  have  chosen  and  not  to  a  limitation  imposed  by  the 
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mathematical  model.  In  order  to  include  the  transient  effect  of  the 


submerged  edge  of  the  floe,  the  finite  difference  grid  should  extend 
several  incident  wave  lengths  in  either  direction  from  the  edge  of 
the  floe.  Further,  to  decay  the  effect  of  the  transient  in  the  verti¬ 
cal  direction,  the  grid  should  extend  several  incident  wave  lengths 
vertically  downward  from  the  edge  of  the  floe.  The  decay  of  the 
transient  in  all  cases  may  be  shown  to  be  of  such  a  form  that  the 
longer  the  incident  wave  length,  the  larger  must  be  the  region 
investigated  around  the  submerged  edge  of  the  floe.  Conversely, 
as  the  incident  wave  length  becomes  shorter,  this  region  diminishes 
in  extent.  Since  the  size  of  this  region  greatly  influences  the 
number  of  grid  nodal  points  investigated  and  hence  the  number 
of  simultaneous  field  and  boundary  equations  (i.e.  the  size  of  the 
solution  matrix),  it  is  clear  that  the  solution  matrix  is  of  a  minimum 
size  for  a  short  wave  length. 

For  the  particular  incident  wave  length  chosen  for  the  numeri¬ 
cal  investigation  conducted  in  this  report  (L  =  10  ft.  ),  we  obtained 
a  solution  matrix  corresponding  to  a  total  of  336  simultaneous 
equations  and  unknowns  which  correspond  to  166  grid  nodal  points. 
The  particular  type  of  matrix  involved  in  this  problem  is  a  banded 
matrix  with  a  band  width  of  about  80.  Such  a  matrix  can  be  solved 
by  the  IBM  7090  using  core  or  immediate  access  storage.  However, 
this  represents  the  practical  limit  of  the  immediate  access  storage. 
Hence,  an  increase  in  matrix  size  or  equivalently,  an  increase  in 
input  wave  length  would  exceed  this  storage  which  would  then 
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necessitate  the  use  of  tape  units.  Although  such  a  computational 
procedure  is  feasible,  the  computer  run  time  could  easily  increase 
by  an  order  of  magnitude  or  more  from  the  run  time  using  only 
immediate  access  core  storage.  Unfortunately,  we  were  limited 
in  this  investigation  to  a  small  computer  budget.  Hence  the  problem 
investigated  was  necessarily  tailored  so  the  resulting  matrix 
solution  did  not  exceed  core  storage.  In  the  event  that  the  com¬ 
puter  run  time  could  economically  cease  to  be  a  governing  criterion 
then  longer  incident  wave  lengths  and  other  ice  floe  thicknesses 
could  be  investigated.  We  feel  that  this  would  be  valuable  from 
several  standpoints. 

First,  the  influence  of  the  incident  wave  length  and  ice  floe 
thickness  on  the  ice  floe  response  could  better  be  understood,  both 
from  a  fundamental  standpoint  and  in  comparison  to  the  neglect  of 
the  finite  floe  submergence.  Secondly,  and  of  equal  importance 
mathematically,  certain  computational  assumptions  used  in  the 
present  analysis  could  be  more  fully  investigated.  These  computa¬ 
tional  assumptions  are  outlined  as  follows: 

1.  In  transforming  the  Laplacian  field  equations  to  finite  dif¬ 
ference  form,  second  order  central  difference  formulas  were 
developed.  However,  along  the  ice-fluid  interface  it  was 
necessary  to  switch  to  higher  order  forward  differences  for 
the  fourth  order  boundary  condition.  This  technique  should 
be  verified  by  increasing  the  order  of  the  finite  difference 
formula  used  in  the  field  equations  to  the  same  order  as  the 


boundary  condition  equations.  Such  an  approach  would  in¬ 
crease  the  width  of  the  banding  in  the  solution  matrix  and 
hence  necessitate  the  use  of  tape  transfer  in  performing  the 
computer  computation. 

2.  There  was  no  check  run  on  the  validity  of  the  choice  of  grid 
size.  Hence  a  computational  run  should  be  performed  on 
the  basis  of  a  smaller  size  grid.  Again,  such  a  process 
would  require  the  use  of  tape  transfer. 

3.  Finally,  the  effect  of  the  region  size  investigated  and  the 
accuracy  with  which  the  computer  satisfies  the  assemblage 
of  simultaneous  equations  should  be  ascertained. 

The  first  two  points,  as  well  as  the  effect  of  region  size, 
can  most  simply  be  answered  by  performing  multiple  solutions; 
varying  certain  parameters  with  each  solution,  and  comparing 
results  in  each  case.  For  instance,  in  the  case  of  the  effect  of  field 
size,  all  other  parameters  being  fixed,  the  field  could  be  increased 
by  varying  amounts  and  the  solution  for  each  case  compared.  When 
an  optimum  field  size  is  reached,  the  solution  at  any  point  within  the 
region  investigated  should  not  further  vary  as  the  field  size  is 
increased. 

The  accuracy  of  computer  solution  can  only  be  checked  by  sub¬ 
stituting  the  answers  back  into  the  original  equations  and  investi¬ 
gating  the  significance  of  the  error  between  the  computed  equation 
value  and  the  input  equation  value.  Such  a  procedure  involves  the 
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analysis  of  solution  stability  and  may  be  investigated  through  the 
use  of  known  techniques. 

However,  the  question  raised  as  to  the  accuracy  of  solution 
will  only  effect  the  quantative  results  of  the  present  analysis.  We 
are,  however,  in  a  position  to  discuss  qualitatively  the  effect  of 
finite  submergence  of  the  floe  for  the  numerical  example  included 
in  this  report.  A  complete  knowledge  of  this  effect  must  naturally 
await  a  detailed  analysis  of  the  influence  of  input  wave  length  and 
floe  thickness  and  an  experimental  substantiation  of  the  mathe¬ 
matical  model  used. 

The  significant  qualitative  results  of  the  numerical  example 
worked,  especially  in  comparison  to  the  results  obtained  from  the 
zero  submergence  case,  are  listed  below: 

1.  Comparing  the  stress  amplitude  variation  in  Figures  4  and 
8  (zero  submergence  and  finite  submergence  respectively), 
it  may  be  observed  that  the  transient  stress  is  about  an  order 
of  magnitude  higher  near  the  leading  edge  of  the  floe  when 
finite  submergence  is  included.  This  result  is  not  surprising 
when  it  is  considered  that  the  barrier  caused  by  the  submerged 
front  of  the  floe  severely  alters  the  local  flow  pattern,  which, 
in  turn,  highly  perturbs  the  local  velocity  characteristics 
inducing  a  pressure  peak  close  to  the  floe’s  leading  edge. 
However,  the  barrier  also  tends  to  greatly  reflect  the  incident 
wave  energy  which  leaves  less  energy  to  be  transmitted  in  the 
form  of  progressive  waves  beneath  the  floe.  Indeed,  a 
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comparison  of  Figures  4  and  8  trnds  to  substantiate  this  in 
that  the  "steady  state"  stress  appears  to  be  of  a  lower  mag¬ 
nitude  in  the  case  of  finite  submergence  than  when  the 
submergence  is  neglected. 

2.  The  transition  to  the  "steady  state"  response  of  the  floe 
appears  to  be  much  more  rapid  in  the  case  of  finite  sub¬ 
mergence.  In  fact,  the  transient  condition  appears  to  die 
out  within  several  wave  lengths  of  the  incident  wave,  indi¬ 
cating  that  wave  length  conversion  and  energy  reflection  is 
taking  place  very  near  the  edge  of  the  floe. 

The  above  observations  lead  one  to  conclude  that  breakup  of 
the  floe  might  well  be  a  progressive  process;  initial  edge  fractures 
which  create  new  boundaries  which,  in  turn,  lead  to  new  edge 
fractures,  etc.  Moreover,  and  of  prime  importance,  is  the  indi¬ 
cation  that  the  magnitude  of  the  steady  state  stress  is  less  in  the 
case  of  finite  submergence.  This  implies  the  possibility  that  a 
conservative  estimate  of  the  bending  stresses  in  the  floe  may  be 
obtained  on  the  basis  of  the  relatively  simple  zero  submergence 
approximation.  The  justification  of  this  must  naturally  await  more 
refined  and  complete  analysis. 

In  NESCO's  opinion,  it  is  important  to  further  pursue  the 
mathematical  investigation  in  the  directions  stated  above.  However, 
it  is  also  very  important  to  obtain  direct  experimental  verification 
of  the  mathematical  model.  Hence,  due  to  the  high  computational 
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costs  attendant  to  a  rigorous  mathematical  extension,  we  feel  that 
the  next  logical  step  in  studying  the  response  of  ice  floes  to  wave 
excitation  is  to  obtain  direct  experimental  evidence  via  field  obser¬ 
vations.  In  particular,  we  urge  that  studies  be  made  of  ice  floe 
deflections  as  a  function  of  time  and  location  in  the  floe.  A  properly 
designed  study  of  this  sort  will  yield  the  period  of  the  waves  present 
in  the  area,  the  magnitude  of  stresses  involved,  their  spatial 
variation,  and  the  dynamic  strength  characteristics  of  the  ice. 
Subsidiary  studies  should  be  conducted  on  the  importance  of  floe 
thickness  variations  and  the  range  of  variation  of  incident  wave 
lengths  and  wave  amplitudes.  A  careful  analysis  of  such  studies 
will  then  indicate  the  proper  direction  to  proceed  mathematically 
and  hence  justify  and  minimize  the  expenses  involved. 

A  detailed  proposal  of  such  an  experimental  study  is  presently 
being  prepared  by  NESCC.  This  proposal  will  be  submitted  in  the 


near  future. 
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